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We present theoretical models, in one and two space dimensions, that exhibit Mott insulating 
ground states at fractional occupations without any symmetry breaking. The Hamiltonians of these 
models are non-local in configuration space, but local in phase space. 



In recent years the study of doped Mott insulators [1] 
has become one of the main themes of condensed mat- 
ter physics [2]. In a Mott insulator the motion of each 
particle is hindered by the presence of other particles 
due to their strong repulsive interaction, and in this re- 
spect they arc different from band insulators where the 
Pauli exclusion principle forbids the motion of electrons. 
An interesting feature of Mott insulators is that doping, 
i.e., changing of particle density, can dramatically affect 
their properties. For example, doping of "cuprates" and 
"mangnites" has lead to the high transition temperature 
superconductors and the colossal magneto-resistive ma- 
terials. 

If a crystalline solid has a fractional occupation num- 
ber [3] and is Mott insulating, there seems always to be a 
spontaneous symmetry breaking of the translation sym- 
metry [4]. After the symmetry breaking the unit cell 
is effectively enlarged so that the occupation number 
becomes an integer. This kind of symmetry breaking 
blurs the distinction between Mott insulators and band 
insulators [5]. 

Nearly thirty years ago, Anderson and Fazekas [6] pos- 
tulated that due to the frustrated antiferromagnetic in- 
teration, the Mott state on triangular lattice does not 
break the spin rotation symmetry (hence is "feature- 
less"). This novel Mott state has since been dubbed the 
name "spin liquid" . Soon after the discovery of high tem- 
perature superconductivity, Anderson proposed that the 
parent (Mott insulating) compounds of these supercon- 
ductors are spin liquids [7]. Although subsequent devel- 
opments have shown this is not true, nonetheless this 
proposal has stimulated wide interest in spin liquids [8]. 

One of the reasons that we are interested in "feature- 
less" Mott insulators is because it is widely believed that 
the elementary excitations in such systems may carry 
fractional quantum numbers (e.g., fractional charge, frac- 
tional spin, etc). As a result, the conducting state(s) 
produced by doping them can be vrey different from the 
those produced by doping ordinary symmetry-breaking 
Mott insulators. 

In this paper we present two examples that exhibit 
featureless Mott insulating ground states. In both cases 
the elementary excitations carry fractional quantum 
numbers. From these examples we learn a valuable 
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FIG. 1: The effect of Hm is to hop a pair of particles while 
conserving their center of mass position. The hopping range 
is set by re -1 . In this figure the nearest neighbor pair hopping 
is illustrated. 

lesson in constructing models of this kind. 

Without further delay let us define these models. First 
we focus on a one-dimensional (ID) model. Consider a 
finite lattice with M = 3N — 2 sites where N is an odd in- 
teger. We label each site by an integer n with the central 
site n = 0. On this lattice we put N spinless fermions. 
The occupation number is given by v — N / (3N — 2) 
which approaches 1/3 in the N — > oo limit. It turns out 
that this choice of the number of sites ensures that the 
ground state is non-degenerate for all N. The Hamil- 
tonian that describes the interaction between these N 
fermions is given by 

H 1D = g K 3 Y,b+b p . (1) 

p 

In the above g is an energy parameter, and is positive. 
The parameter re, whose meaning shall be discussed lat- 
ter, is dimensionless (0 < re < 1). Since Hid is a sum 
of positive-definite operators, its eigenvalues are clearly 
non-negative. 

In Hid the index p runs through both integers and 
half integers, and 

K = ^2{ie~ K2q2 )c p - q c p+q , (2) 
i 

where q takes integer/half- integer value depending on 
whether p is integer/half-intcgcr. In Eq. (2) c p ± q an- 
nihilates a fermion at lattice site p ± q. Hid describes 
the hopping of pairs of fermions from sites p — q and p + q 
to p — k and p + k (Fig.l). In such a hopping process the 
center-of-mass position (jp) of the fermion pair remains 
unchanged. From Eq. (2) it is apparent that re -1 con- 
trols the hopping range. 

For occupation v = N/ (3N — 2) and any re, the model 
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defined by Eq. (1) possesses a ground state given by 
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This state is annihilated by all the bp's, i.e. 

bp^iD > = 0. 



(4) 
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As a the result Hid\^id >= implying \^id > is a 
ground state of H\d- Eq. (5) can be shown to be true 
by rewriting it in terms of the wavefunction x({ n j}) an( i 
using the following identitiy[9] 
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The $ p in Eq. (6) is given by 
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In latter part of the paper we shall show that for suf- 
ficiently small positive values of re, \^Sid > is the unique 
ground state of Hid , it separated by an energy gap from 
the lowest excited states, and it does not break transla- 
tional symmetry. Meanwhile, to get a feeling for Y&id >, 
let us look at a representation of it for 3 particles, as 
shown in Fig. 2. The result is a coherent superposition of 
many different configurations. The relative weight and 
phase of these configurations ensures that they are anni- 
hilated by bp. Any breaking of the weight-phase relation 
causes excited states. 

Next we present a two dimensional (2D) model with 
similar properties. The lattice consists of M' chains, with 
each chain made up of M = 3N — 2 sites. In each chain 
we place N spinless fermions. As the result, the 2D oc- 
cupation number is also given by v — N/(3N — 2). In 
the limit N, M' — > oo the occupation number is also 1/3. 
The Hamiltonian of the 2D model is given by 
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FIG. 2: \ipiD > for three particles. It is represented by 
components with different occupation of sites. The relative 
weight and phase between these configuration are e 9ft (top), 
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In the above the chain index n runs through integers and 
p can be both integer and half- integers. In H 2 r> 
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where q runs through integer/half-integer depending on 
whether p is integer/half-integcr. The energy parameters 
gi and g 2 are both positive. The first term of H 2 d hops a 
pair of particles within each chain n, and the second term 
hops a pair of particles from chain n to chain n + 1 and 
vice versa. In both hopping processes the center of mass 
coordinate along the chain direction remains unchanged. 
For any values of g\ and g 2 , 
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where \^\D, n > is the previous \^id > defined for the 
nth chain, is a zero-energy eigen state of H 2 £>. This is 
because b p annihilates \^id >■ Let us imagine obtaining 
H 2 d by turning on g 2 from zero while holding g\ fixed. 
At non-zero g 2 the excited states of adjacent chains are 
mixed. However, so long as g 2 is small compared to the 
excitation gap at 172 = 0, this mixing will not produce 
negative energy eigen states. [10] Under that condition 
\1p2D > continues to be the exact ground state of H 2 jy, 
and it is separated by an energy gap from the excited 
states. Since \ipiD > does not break translation symme- 
try, so does'nt \ 1 $' 2 d >• 

The remaining task is to collaborate the statements 
we have made about \ipiD >■ To accomplish this we 
will make use of a mapping from the ID Hamiltonian to 
a well-studied problem in two dimensions, namely, the 
spinless fermions in a strong magnetic field. 

It is well known that when the lowest Landau level 
(LLL) is 1/3 filled, and if the interaction between 
fermions falls off sufficiently fast as a function of the in- 
terparticle distance, the ground state is an incompressible 
quantum fluid. It is also well-known that when a particle 
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is confined to the lowest Landau level the two compo- 
nents of its coordinate no longer commute. This implies 
that the dynamics of the fermions is effectively that of 
a ID system. This intrinsic one-dimensionality can be 
made explicit through a mapping of the LLL problem to 
one dimension. What has not been appreciated is that 
although the time reversal symmetry is explicitly broken 
in the 2D problem, the corresponding ID Hamiltonian is 
time-reversal invariant. Moreover, after the mapping to 
ID, the ground state is an example of a featureless Mott 
insulator at 1/3 occupation number. 

To start, let us consider a spin-polarized two dimen- 
sional electron gas where the electrons interact through 
a two-body potential 

H int = V J d 2 rd 2 r'V(r-r')ij + (r)^ + {r')^(r')^(r) 
V(r) = V 2 <5(r). (11) 

Although this interaction potential is a singular function 
of the coordinates, it has non-singular matrix element 
between states in the lowest Landau level. Laughlin's v — 
1/3 wavefucntion [11] is known to be the exact ground 
state of this potential [12]. The gap to the lowest excited 
state is proportional to the parameter Vo/l%, where Ib is 
the magnetic length. 

Next we place the above two dimensional electron gas 
on a cylinder with circumference L with the magnetic 
field perpendicular to the surface. In the Landau gauge 
the LLL basis orbitals are of the following form 
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These orbitals are delocalized around the cylinder (y) 
while localized along its axis (x). The quantization of 
the y- momentum in units of 2tt/L explains why we have 
a lattice. If we expand the field operator tp(r), projected 
to the lowest Landau level, as V( r ) = £n^«( r ) c ™' su b~ 
stitute the result into Eq. (11), and rescale the coor- 
dinates by L/2w, we obtain the Hamiltonian (1) with 
k = 2-kIb/L, and g = w^jiji " Thus k measures the 
ratio between the magnetic length and the circumference 
of the cylinder. Since k^ 1 controls the hopping lengths 
of the ID and 2D Hamiltonians , it is clearly desirable 
to restrict k > 0. On the other hand, to preserve the 
the liquid property of the Laughlin state we need to keep 
L > Ib (or k < 2ir). To satisfy both requirements we 
assume < k < 1. 

In Landau gauge, the v = 1/3 Laughlin wavefunction 
on a cylinder reads [13, 14] 
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where z = x + iy. The x({ n k}) given in Eq. (4) are the 
coefficients of when expanded in terms of products 
of the lowest Landau level orbitals (Eq. (12)) . 



At this point we have mapped Hi nt to H\d, and the 
Laughlin state |^i/3 > to \$>id >• Since the above map- 
ping is a unitary transformation, all known properties of 
1^1/3 > and H int are preserved. For example, the fact 
that I ^1/3 > describes a quantum liquid at magnetic fill- 
ing factor 1/3 translates into the statement that \^\d > 
is a quantum liquid at occupation number 1/3. The fact 
that |*i/3 > is the non-degenerate ground state of Hi nt 
implies that [^id > is the non-degenerate ground state of 
Hid- Finally, the fact that Hint possesses an energy gap 
and fractional charge quasiparticles at filling factor 1/3 
implies the same for Hid at lattice occupation 1/3. In 
this way we have established that the model presented 
by Hid possesses a featureless Mott insulating ground 
state at fractional occupation number! 

By using the ID to 2D mapping discussed previously 
it is simple to map H2D to a 3D Hamiltonian H' int de- 
scribing coupled layers: 



W 

n int 



n 

H n = yd 2 rd 2 r'F(r-r')[^^+(r)V+(r')^(r')^n(r) 
+ ^ 2 (^+ +1 (r)V+ +1 (r')^(r')^(r) + /i.c.)]. (14) 

In the above equation V(r) is the same as that given in 
Eq. (11), and oc g\^ respectively. With V2 = the 
problem becomes that of many independent layers. At 
v = 1/3 each layer is in the liquid ground state described 
by I* 1/3 >■ The term proportional to V2 hops a pair 
of electrons between adjacent layers. As discussed below 
Eq. (10), so long as V2 is sufficiently smaller than V\, the 
ground state is the direct product of the 1/3 Laughlin 
liquid in each layer. 

It is interesting to note that unlike usual Hamiltonians, 
Hid and H2D do not have the single-fermion hopping 
term. One can in fact add a single particle hopping term 



H' 1D = -te- 2 /^ [4 +1 c n + h.c] 



(15) 



to Hid without changing qualitatively any of the results. 
To understand this we note that after mapping to 2D 
Eq. (15) becomes 



-t J d 2 rcos(2ny/L)tp + 



(r)^(r), 



(16) 



i.e., a periodic potential. It is clear that so long as |t| 
is much smaller than the energy gap of Hi nt , it will not 
change the quantum liquid nature of the ground state 
and will not collapse the excitation gap. Similar single- 
particle terms can by added to H 2 d- Again, so long as 
the strength of the single particle hopping terms is suffi- 
ciently small, the qualitative results we discussed above 
will be preserved. 
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Models that gives featureless Mott insulating ground 
state at fractional occupation number have been difficult 
to find. What is special about the models presented in 
this paper? To answer this question let us focus on Hid . 
A standard (ID) interaction potential is local in config- 
uration space (q-space), i.e., 

V = YlJ dq ^U(q t - qj)\qi, ..,q N >< qi, -,?iv| (17) 

3 i<3 

When this type of potential dominates the Hamiltonian, 
it favors particles to assume a fixed q-space configuration. 
For a repulsive potential such a configuration often take 
the form of a regular lattice, ff^, on the other hand, is 
non-local in q-space and the particles do not necessarily 
have to pay a high price in energy to get close. However, 
there is a hidden locality in Hid which is revealed not in 
configuration space, but in phase space. 

The quantum description of phase space dynamics is 
most conveniently done in terms of coherent states. For 
a single particle with coordinate q and momentum p, the 
coherent state \z > satisfies 



X ■ 

q + i—p 



z >-- 



z \z >, 



(18) 



where z = x + iy and A is any pre-chosen length scale 
serving to make q and (X 2 /h)p the same dimension. If 
\q > denotes the position eigenstate, it is simple to show 
that 



< q\z > oc e 
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and this function defines the transformation between the 
phase space and configuration space descriptions. It is 
interesting to note if we identify A with the magnetic 
length and q/\ 2 with the y- momentum, the functions 
given in Eq. (19) become precisely the LLL basis orbitals 
in Eq. (12). Thus a 2D quantum mechanical problem de- 
fined in the LLL is completely equivalent to the coherent 
state representation of a ID quantum problem. 

Let us go back to the non-local Hamiltonian Hid- 
When expressed in terms of the coherent state basis, it 
gets a local form 



to phase-space local potential is the main reason that 
they have featureless Mott insulating ground states at 
fractional occupation numbers. 

In conclusion, we have constructed ID and 2D lattice 
models that exhibit featureless Mott insulating ground 
state at fractional occupation number. The Hamiltoni- 
ans of these models have finite range interactions, and 
respect lattice translation and time reversal symmetries. 
These Hamiltonians are local in phase space instead of 
configuration space. We believe phase space local inter- 
action could be a missing key for constructing models of 
this type of novel insulators. 
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Zj)\zi, .., z N >< zi, .., z N [20) 



The V(zi — Zj) in the above equation is the potential 
given in Hi nt . This potential prevents particles from get- 
ting close together in the phase space. However due to 
the non-commutativity of the phase space coordinates, 
the particles are not frozen in any particular configura- 
tion. We believe the fact that Hid and H 2 d are linked 
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